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Abstract. A free boundary problem describing small deformations in a membrane based model of electro- 
statically actuated MEMS is investigated. The existence of stationary solutions is established for small voltage 
values. A justification of the widely studied narrow-gap model is given by showing that steady state solutions 
of the free boundary problem converge toward stationary solutions of the nan'ow-gap model when the aspect 
ratio of the device tends to zero. 



1. Introduction 

Microelectromechanical systems (MEMS) have become key components of many commercial systems, 
including accelerometers for airbag deployment in automobiles, ink jet printer heads, optical switches, 
micropumps, chemical sensors and many others. Idealized modern MEMS devices often consist of two 
components; a rigid ground plate and a thin and deformable elastic membrane that is held fixed along its 
boundary above the rigid plate, and its design is based on the interaction of electrostatic and elastic forces. 
More precisely, when a voltage difference is applied between the two components, a Coulomb force is in- 
duced which is varied in strength by varying the applied voltage and gives rise to deformations of the elastic 
membrane. Perhaps the most ubiquitous nonlinear phenomenon associated with electrostatically actuated 
MEMS devices is the so-called "pull-in" instability limiting the effectiveness of such devices. In this insta- 
bility, when voltages are appUed beyond a certain critical pull-in voltage, there is no longer a steady-state 
configuration of the device where the two components remain separate. This possible touchdown of the 
membrane on the ground plate affects the design of the devices as it severely restricts the range of stable 
operation. The understanding and control of the pull-in voltage instability are thus of great technological 
importance: in this connection, a large number of MEMS devices which rely on electrostatic actuation have 
been investigated both experimentally and through numerical simulations and several mathematical models 
describing these devices have been set up. 

We consider here a simple membrane based model of an electrostatically actuated MEMS device as de- 
picted in Figure [T] and refer the reader e.g. to |fT9l l20l l22l and the references therein for a more detailed 
account of the physical background and the modeUng aspects of modern MEMS devices. In this simpUfied 
situation, we assume that the applied voltage and the permittivity of the membrane are constant (normalized 
to one) and that there is no variation in the horizontal direction orthogonal to the a--direction of both the 
(dimensionless) electrostatic potential tp and the displacement u of the membrane. Under appropriate scal- 
ings, the rigid ground plate is at z = — 1 and the undeflected membrane at z = is fixed at the boundary 
a; = — 1 and x = 1, see Figure [T] Denoting the aspect ratio of the device, i.e. the ratio of the undeformed 
gap size to the device length, before scaling by e, the dimensionless electrostatic potential ip — ip{x, z) is 
supposed to satisfy Laplace's equation 

e^dl^ + dl^ = (l.I) 

in the region 

niu) := {(x, z) e (-1, 1) X (-1, oo) : -1 < z < u{x)} 
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Figure 1 . idealized electrostatic MEMS device. 



between the rigid ground plate at z = — 1 and the deflected membrane at z = u. The boundary conditions 
are then 

^ ^ on z -1 (1.2) 

and 

ip = 1 on z = u . (1.3) 

As for the deformation of the membrane, it results from a balance between dynamic, electrostatic, and 
elastic forces and the membrane displacement u = u{t, x) £ (—1, oo) evolves according to Newton's law 

a^d^u + dtu-dlu = -X {e"^ \da,i;{x,u)\^ + \d^ij{x,u)\'^) (1.4) 

with clamped boundary conditions 

it = at x = ±l. (1.5) 

In ( 11.4b . the term dtu and the right-hand side account for a damping force and the electrostatic force, 
respectively, while the term d^u describes the deformation due to stretching. The latter is obtained after 
linearization resulting from the assumption of small deformations. The contribution to deformation due to 
bending may also be included in ilAi by adding a fourth-order term Bd^u, B > 0, to the left-hand side 
of ( 11.4b but is neglected here. The parameter A > characterizes the relative strengths of electrostatic and 
mechanical forces. It acts as a control parameter proportional to the applied voltage. The coefficient a > 
is indirectly proportional to the damping coefficient. 

Observe that (11. lb is a free boundary problem as the domain between the rigid ground plate and the 
elastic membrane changes with time. Due to this, equations dl.lb and (1 1.4b are strongly coupled. However, 
a common assumption made in all mathematical analysis hitherto is a small aspect ratio e. Formally, sending 
e to zero allows one to solve explicitly ( ll.lb -( ll.3l ) for the potential ip — ipQ, i.e. 

I + z 

V'o(a;, z) = -— — — , {x, z) € Q.{u) , (1.6) 
1 + u(x) 
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thus reducing the free boundary problem to the small aspect ratio model, that is, to an evolution equation 

a^dfu + dtu-dlu = ~—^^, e (0,oo) X (-1,1), (1.7) 

(1 + u)^ 

subject to ( 11.51 ) solely involving the displacement u. Note that the boundary conditions il.2i . il.3i suffice 
to determine which then satisfies on the lateral boundaries a; = ±1 

^/'o(±l,z) = 1 + z, ze(-l,0), (1.8) 

due to il.5\ . The small aspect ratio model with ( ll.SI l has widely been investigated in the recent past 
(with possibly an additional fourth-order term accounting for the deformation due to bending as already 
discussed) and also variants thereof, e.g. in higher dimensions or with additional permittivity profiles or 
non-local terms. An obvious difficulty arising in the study of (11. 7t is the singularity of the source term 
— A/(l + it)^ as u approaches —1 which corresponds to the aforementioned touchdown phenomenon for 
the MEMS device. Concerning the dynamic behavior of small aspect ratio models we refer the reader to 
E [13] [HI for the hyperbolic case a > and to E H] [TT] [HI [SI [ID for the corresponding parabolic 
equation with q = when damping or viscous forces dominate over inertial forces in the system. The 
dynamic behavior of a membrane evolving according to ( 11.7b . (11.51 1 is determined by a pull-in voltage 
A* > 0, see e.g. 13] [7] [TT] [ID . More precisely, if A < A, there are a stable and an unstable steady state (i.e. 
time independent) solution of ( 11.7b subject to (11.5b . that is, of 

dlu ^ ^ xe(-l,l), ui±l) = 0, (1.9) 

(1 + u)^ 



and solutions to the dynamical problem (11.7b starting out from u = converge toward the stable steady 
state. Steady states cease to exist for voltage values A above A* and solutions to the dynamic problem touch 
down on the ground plate in finite time, that is, a pull-in instability occurs. We refer the reader to B1[5l l20l 
for a review of these results and references as well as to |[3]|6]|7]|8][TT|[12][T9]|2T| and the references therein 
for further details on small aspect ratio models. 

To the best of our knowledge, the original model without small gap assumption has not been tack- 
led so far from an analytical point of view. The aim of this paper is to make a step in this direction 
by studying the stationary free boundary problem. More precisely, we shall focus on finding functions 
u : [—1, 1] —J- (—1, oo) and ip : n{u) — > R satisfying the coupled system of elliptic equations 

s'^dltlj{x,z) +d^,^p{x,z) =0 , {x,z)€n{u), (1.10) 

tpix^z) = 1 + z , {x,z)edn{u), (1.11) 

dlu{x) = X{e^\d,,^{x,u{x))\'' + \dMx,u{x))\^) , xG(-l,l), (1.12) 

u{x)^0, x = ±l, (1.13) 



where the domain of definition n{u) of ijj is 

Q{u) {(x, z) ; -1 < X < 1 , -1 < z < u{x)} , 

that is, the two-dimensional region between the rigid ground plate and the membrane with deflection u. 
Obviously, f2(u) is a domain and possesses four corners provided the values of the continuous and convex 
(see (11.12b ) function u satisfying ( 11.13b stay away from —1. Let then 

r{u) dn{u) \ {(±1, -1), (±1, 0)} 

denote the boundary of il{u) without corners. System (I1.10b - (ll.l3b is exactly the time-independent version 
of ( |l.lb - (ll.5b subject to the lateral boundary condition dl.Sb which is imposed to make the system well- 
posed. This particular choice of a continuous boundary condition is made mainly for simplicity and we 
point out again that this condition is satisfied by ipQ from ( |1.6b in the small aspect ratio limit. For the 
stationary free boundary problem we shall show existence of smooth solutions for small voltage values A: 
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Theorem 1.1. There exists Ag > independent of £ € (0,1) such that (|1.10t -( fI.13t admits for each 
A G (0, Aq] a solution 



u e c2+"([-i, 1]) , V e w|(fi(w)) n c{niu)) n c^+"{n{u) u t{u)) , 

where a G [0, 1) w arbitrary. The function u is even, convex, and satisfies 

0>u{x)>-1 + kq, a-G(-l,l), (1.14) 

\\u\\wu-i-i) < — ' (1-15) 

Kq 

for some kq G (0, 1) independent of e. Moreover, ip = ipix, z) is even with respect to x G (—1, 1)- 

We refer to Section|2]for the proof of Theorem I 1 . 1 1 which is based on a transformation to a fixed domain 
and on an appUcation of Schauder's fixed point theorem. Concerning the latter, given a displacement 
u G W^{—1, 1) with values in (—1, 0), we first construct in Lemma l272l the corresponding solution t/ju to 
dl.lOb - dl.lll l by using an equivalent formulation on a rectangle. Of particular importance is the regularity 
of the trace of the gradient of t/ju on the upper boundary {z = u} as stated in Lemma 1241 which plays an 
important role in the subsequent analysis of ( ll.l2b -( fTTT3] l. Indeed, it is used as a source term to construct 
a solution S{u) to ( ll.l2l )-( fl~T3] ) with instead of tp, see Lemma [23] Restricting suitably the set of 
admissible displacements u, the map S turns out to enjoy the properties needed to apply Schauder's fixed 
point theorem. 

In particular, for values A < Aq, Theorem I 1 . 1 [ provides for each e G (0, 1) a solution {u^, to dl.lOl )- 
( |1.13t satisfying the bounds ( |1.14t . ( 11.151 ) uniformly with respect to e G (0, 1). This property allows us 
to give a rigorous justification of the small aspect ratio model (11. 7b . (11. 5b by showing that (ug, V'e)ee(o,i) 
converges toward a solution to that model as e tends to zero. More generally, we have: 

Theorem 1.2. Let A > and let {u^,'4'e)e<^(o,i) be a family of solutions to (I1.10b -( ll.l3b satisfying the 
bounds (11.14b and ( 11.15b . Then there are a sequence \ and a (smooth) solution uq to the time- 
independent small aspect ratio equation (11. 9b such that 

— >uo in M^i(-l,l) 

and 

^EfclflKJ ^V'olj2(uo) i2((-l,l) X (0, 1)) (1.16) 

as k oo, where V'o th^ corresponding potential ( 11.6b with u = uq. 

The proof of Theorem 1 1.21 is performed in Section |3] by using a compactness argument. Since ( 11.101 ) 
becomes degenerate elliptic in the limit e — > 0, the regularity of ipe is no longer the same in the x- and 
^-directions and a cornerstone of the proof is to obtain estimates for the trace of dzipe on {z = u^}. 

2. Existence for small voltage values: Proof of Theorem 1 1.1 1 

We first prove Theorem ll.il Since the domain of definition of the potential t/' in (II. 10b depends on the 
displacement u of the membrane, we use an alternative formulation by transforming the problem on a fixed 
domain, that is, on the rectangle O := (—1, 1) x (0, 1). More precisely, given a function u G W^{—1, 1) 
taking values in (—1, oo) and satisfying the boundary conditions u(±l) = 0, we define a diffeomorphism 
Tu := ^l{u) Qhy setting 

l + z 



Tu{x,z) :^ , {x,z)&n{u). (2.1) 

V l + u{x)J 

Clearly, 

T-\x,Tj)^{x,{l + u{x))i^-l) , (x,77)Gn, (2.2) 
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and it readily follows that problem ( ll.lOt - dl.llt is equivalent to 

(x, 77) = , {x,T]) en , 

(j){x,rii) = 77 , (a;,?/) G dVl , 
for = 1/) o T^^, the u-dependent differential operator being defined by 



(2.3) 
(2.4) 



CuW :=£ d^w — 2e''7j 

dxu{x) 



2 02 „ 



1 + u{x) J 1 + u{x) 



drfW 



Moreover, ( |1.12t . dl.lSt become 



dlu{x) = A 



(-1,1), 

x = ±l, 



(2.5) 
(2.6) 



{l + u{x)Y 

u{x) = , 
where we have used 

1) , xe(-l,l), (2.7) 

since 0(.t, 1) = 1 for x e (-1, 1) by ( 12.41 ). 

Our goal is to solve ( I2.3l )-( l2.6l l by means of Schauder's fixed point theorem. Fixing ro € (0, 2), we 
introduce the set 

C := {ueW^{-l,l)nWli){-l,l) : u is even and < < tq} , 

where 

Wljji-1,1) {ueW^i^l,!) ■■ u{±l) = 0} , [1,00]. 

Let us first collect some properties of C. 

Lemma 2.1. C is a closed, convex, and bounded subset 0/ 1, V) for each g € [1, 00] and 



> > -y > -1 , a-e(-l,l), ueC 



(2.8) 



Proof. Clearly, C is convex and closed in 1, 1) and thus weakly closed in W^(— 1, 1). Therefore, C 

is convex and closed in Wg{~l, 1) for each q G [1, 00]. Next, for u e C, integrating the equality 

dxu{x) = dxu{y) + [ dlu{z)dz, (a;, j/) £ (-1, 1) x (-1, 1) , 



with respect to y on (—1, 1), we find: 

\dM^)\<2ro, xei-1,1). (2.9) 

Since u(±l) = 0, we deduce from ( 12.9b that \u{x)\ < 2ro for x E [—1, 1] and C is thus bounded in 
W^l—l, 1). Next, the convexity of u and the boundary values u{±l) — clearly ensure that u < 0. 
Finally, if u attains a negative minimum at some point Xm G (—1, 1), we may assume Xm G [0, 1) without 
loss of generality since u is even. Then dxu{xm) = and 



u{x)-u{xjn)= / dxu{y)dy = [{y-x)dxu{y)] 
Thus, since u(l) = 0, 



y=x„ 



{y~x)dlu{y)dy 



{x-y)dlu{y)dy . 



-u{Xm) 



{l-y)d^^u{y)dy<^ , 



from which (12.81 1 follows. 

Next we study the existence and properties of the solution to ( I2.3l l- (I2.4| | when w G C is given. 



□ 
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Lemma 2.2. Given u <eC there is a unique solution (/>„ g 14^2 (^) to ( |2.3l )-( |2r4l ). Moreover, 0„ = (/)m(x, 77) 
is even with respect to x. 



and 



ri{l + u{x)) < cj)u{x,ifii) <1 , {x.ifj^Vt, uGC, 



(2.10) 
(2.11) 



for some constant ci = ci(ro,e) > 0. 



Proof. We claim that the operator —Cu is elliptic for u E C given. To see this, choose an arbitrary u E C 
and let 



A 



1 + u{x) 



V 



z^dxu{x) 1 + e'^Tf{dxu{x)Y 



V l + u(.r)'' + / 

denote the principal part of — £„ for fixed (a;, 77) £ f2 with trace t and determinant d given by 



t := e 4 



(1+U(x))2 



Then the two eigenvalues of A are 



/x± = i(i±v/i2-4d) 



and since 



by (El and dZOl ), 



1 + £^ < t < £^ 



/J+ > M- ^ — 5^ 



fl + 4e2r2) 



> 0. 



t '- e2(2-ro)2 + 4+16e2r2 

Consequently, — £„ is elliptic with a positive ellipticity constant depending on ro and e but not on u e C. 
Next observe that (I2.3b -( l24l i is equivalent to 

(£.$)(x,7?) = -/„(a^,ry), (x,77)en, (2.12) 

^{x,r]) = 0, ix,rj)€dn, (2.13) 



by setting 77) := (j){x, rf) ~ 77, (x, 77) £ Vl, where G Loo(f^) is defined as 



fuix,T]) := £„77 = £-^77 



dxu{x) \ d^u^x) 



and satisfies 



11/ 



< e 



1 + ^(a::) / 1 + u{x) 
2ro 



{x, r,)en 



(2.14) 



(2.15) 



,(2-ro)2 2-ro, 

by ( 12. 9t and Lemma IZTl Noticing that, thanks to Lemma IZTl all coefficients of £„, written in divergence 
form 



CuW = dx [ e^dxW - e^rj 



dxu{x) 
1 + u{xy 



d^w + a,, -£^77 



dxu{x) 
'l + u{x) 



dxW 



l + £^ri^{dxu{x)) 



{l + u{x)) 



2 drjW 



2 dxu{x) 2 
£ -T- — r^OxW - e -q 



1 + u(a;) 



dxu{x) 
1 + u{x) 



driW , 
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as well as /„ have norms in Loo{il) uniformly bounded with respect to u e C, we may apply |fT6l Chapt. 3, 
Thm. 9.1 & Thm. 10.1] to obtain the existence and uniqueness of a solution G W2 to (I2.12l l- (l2.13l l 
satisfying 

<c(||$„|U,(o) + l) (2.16) 

with a constant c depending on ro and e but not onu E C. Setting (puix, rf) = $m(.t, rf) + i] for (.t, rf) E Cl, 
the function (pu obviously solves ( I2.3l )-( l2.4b and, owing to ( 12.16b . the bound ( 12.1 lb readily follows provided 
we can verify (12.10b . For this we take w = 1 and note that CuW = in 51 while 1(7(77) = = 'Puix, rj) 

for {x, rj) € d^l. The comparison principle then ensures 0„ < 1 in fj. Taking v{x, rj) :— ri{l + u{x)) for 
{x, rf) € f2, we have CuV = in and v{x, 1]) <ri = 4>u{x, rf) for (x, rf) £ d^l. We conclude that (f>u > v 
in f2 again by the comparison principle, and ( 12.10b and ( 12.1 lb follow. It remains to check that (pu is even. 
However, u being even, dxU is odd and d^u is even, and it is easily seen that (j>{x, rf) := (f>u{—x, rf) satisfies 
( I2.3b -( l24l ) as well, whence (p = by uniqueness. □ 

We next turn to the continuity property of with respect to u G C. 

Lemma 2.3. The mapping (u i— > 0„) : C — > (^) '■^ continuous when C is endowed with the topology 
ofWii-1,1). 

Proof. Let C be endowed with the topology of 1, 1). Given u G C, we define a bounded linear 

operator ^(u) G C{Wljj{n), L2in)) by 

:=-£„$, $GW|,^(f7), 

and note that A is continuous from C in £ (W| ^^(ri), £2(^1)), thanks to the continuous embedding of 
Wii-1, 1) in W^{-1, 1) and the boundedness'of C in W^^l-l, 1). Then CS Chapt. 3, Thm. 9.1 & 
Thm. 10.1] (see the proof of Lemma |272] | warrants that A{u) is invertible for each u E C. Owing to the 
continuity (in fact; analyticity) of the inversion map £ i-> of bounded linear operators, we conclude that 

C ^ C{L2{fl),Wljj{n)) , u^Aiu)-^ 

is continuous. One then checks that u i-> fu is continuous from C to ^2(^1), where /„ is given in ( 12. 14b . 
Consequently, 

(u ^ - A{u)-\-fu)) : C Wlom 
is continuous. Recalling that cf>u{x, rj) = ^u{x, rf) + 77 for (x, rf) E^ gives the claim. □ 

To obtain estimates on solutions to ( ll.l2b -( fTTT3] l we need estimates on the gradient of (j)u on the boundary 
77 = 1 as provided by the following lemma. 

Lemma 2.4. There is a constant C2 > depending only on tq G (0, 2) and e G (0, 1) such that, given 
u £ C, the corresponding solution 0„ G M^2 (^) ( I2.3b -( I2.4I ) satisfies 

\\dM;mwl/-(.,^^)<C2 (2.17) 

and 

< 1) < 1 + 2e^ a;G(-l,l). (2.18) 

Proof. According to jlSl Chapt. 2, Thm. 5.4] there is a positive constant c depending only on Q, such that 

\\dn't'u{-A)\\wl/^{-i,i) ^ c||0„||H^2(n) 
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from which (|2.17| l readily follows by ( 12.111 ). Next, set Waif]) := 77^+" for 77 G [0, 1] and a > 0. Then 

Wa{v) < = fpuix^ 1]) for {x, rj) e dQ and 



(1 + uj^ 



dxU \ ^ d^u 



1+u 1+u 



~ (1 + 1t)2 L ' ^ ^ -I - (1 + u)2 L -I 

in n, where we used u < and < d^u < ro < 2 to obtain the last inequality. Consequently, choosing 
a = 2£^, we realize that £„W2e2 > in fJ, and we infer from the comparison principle that 

In particular, for 77 S (0, 1), 

' K{X,T]) - (f>u{x,l)) = {(t)u{x,r]) -1) < ^—-{W2e^{v) - 'W2e'2{l)) , 



?/ — 1 r/ — 1 

whence dri4>u{x, 1) < driW2e^ (1) =1 + Se^ for a; e ( — 1, 1). Since 0„ < 1 in and (j)u{x, 1) = 1, we also 
havea,,(/)„(x,l) > Ofora; G (-1,1). □ 

Next, given u e C, we set 



i+f{dMx))^ 

{l + u{x)Y 



9u{x) := ,^ ; ; ,\;^ \d^^u{xA)\' , a-e(-l,l), (2.19) 



and observe that Lemma \TA\ and ( 12.181 ) guarantee that G Loo(— 1, 1)- Thus, for each A > there is a 
unique solution v = S'(u) in _d(~1: 1) to the linear problem 

dlv{x) = \gu{x) , a; e (-1, 1) , (2.20) 

w(x)=0, x = ±l. (2.21) 

Actually, we have: 

Lemma 2.5. IfC is endowed with the topology 0/ Wf (— 1, 1), then S : C VFf ~'''^(— 1, 1) is continuous 
for each a £ [0,1 /2), and there is a positive constant C3 (cr) depending only tq, £, and a such that 

II'5'(")IIm/2+,.(_j -^^ < A C3(cr) , mgC. (2.22) 
Moreover, S{u) is even and convex for u G C. 

Proof. Lemma |23] together with ifTSl Chapt. 2, Thm. 5.4] and ( 12.171 ) imply that u ^ dri4>u{-, 1) is contin- 

1/2 

uous and bounded as a mapping C W2 (—1, 1). In the following, given two Banach spaces X and Y 
of real-valued functions, we write X ^ Y to indicate that X is continuously embedded in Y and we set 
X -Y := {fg : (/, g) & X x Y}. Let < a < ai < 1/2. Since pointwise multiplication 

w^/\-i, 1) • w.^/\-i, 1) ^ i^r (-1, 1) 

is bilinear and continuous according to 12] Thm. 4.1 & Rem. 4.2(d)], we infer that u H> l)p defines 

a bounded and continuous mapping C — > (—1,1). Noticing that 

C~>W^{-1,1), u^— — (2.23) 

(1 + u{x))^ 



is continuous and bounded as well by Lemma 12711 and that 

W^{-1, 1) • W^' (-1, 1) ^ W^{-1, 1) 
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by 12] Thm. 4.1 & Rem. 4.2(d)], we see that u i-> g„ is continuous and bounded from C to VF^f (— 1, 1). 
Consequently, S : C — > W2^'^(— 1, 1) is continuous and satisfies (|2.22| i. Clearly, S{u) is even for u <E C 
since u and 0u(-, 1) are even by Lemma lZSl and S{u) is convex since d'^S{u) > by ( |2.19t and ( I2.20l i. □ 

We are now in a position to construct solutions to ( |2.3t -( |276] l for small values of A by applying Schauder's 
fixed point theorem to the map S. 

Proposition 2.6. There exists Aq > independent of e E (0, 1) such that (I2.3b -( l2.6l l admits for each 
A G (0, Ao] a solution 

satisfying 

> u{x) > -y > -1 and < d'^u{x) < vq , x e (-1, 1) . 
Moreover, u is even and belongs to W^^" (-1,1) for any a e [0,1/2). 

Proof. To prove that S maps the closed and convex subset C of 1, 1) into itself note that ( 12.91 ), ( 12.181 ) 

together with Lemma lZTI ensure 

< dlS{u) = X ^+fff[f" \d,M-, 1)P < 4Al±i^(l + 2e') (2.24) 
{l + u{x)y [2-roy 

for u £ C. Thus there is Ao = Ao(ro) > sufficiently small and independent of e G (0, 1) such that 
< dlS{u) < ro for A £ (0, Ao] and u G C, so it follows from Lemma [275] that S indeed maps C into 
itself. Since W2^+'^(-l, 1) embeds compactly in W^{-1, 1) for cr e (0, 1/2), Lemma [23] impUes that 
S : C —i' C is continuous and compact and thus has a fixed point u E C enjoying the properties stated in 
Lemma 12. 11 □ 



Clearly, a positive lower bound on Aq can be obtained by optimizing its choice according to ( 12.241 ). 

To finish off the proof of Theorem ll.ll it remains to improve the regularity of {u, (pu) and to pull it back 
on the domain ri(w) by means of the transformation Tu from (12.1b . 

Corollary 2.7. If (u, 0„) is the solution to (|2.3b -( l276] l for A G (0, Aq] provided by Proposition \2.6\ then 
{u, ip) with ip :^ (pu ° Tu is a solution to ( |1.10b -( [I.13b with regularity 



u e 6*2+" ([-1, 1]) , V e w^{n{u)) n c{niu)) n c^+''{n{u) u t{u)) , 

for each a G [0, 1). 

Proof. From (/)„ G M^2^(ri) and ( 12.11 ) we readily deduce = o G M^2^(r2(u)) and solves dl.lOb 
in H{u). Moreover, since il{u) is a Lipschitz domain, the trace of is well defined as an element of 
W2^^{dn{u)) according to Ql Chapt. 2, Thm. 5.5] and ( [mi l follows from u{±l) = and Also, 
since satisfies the exterior cone condition at every point of its boundary and u G W^{—1, 1), it follows 



from Thm. 9.30] that G C{n). RecalHng that G C {niu); Qj, we deduce that ^/j e C [n{u 
Finally, -0 is even in x due to Lemma [2721 and the fact that u is even. 

We next improve the regularity of ijj with the help of |[T0] Thm. 5.2.7]. To this end, we note that, 
since u G W^{—1, 1), the boundary dil{u) of il{u) is a curvilinear polygon of class C^'^ in the sense of 
|[T0] Definition 1.4.5.1] with four vertices {(-1, -1), (1,-1), (-1,0), (1,0)} connected by W^^-smooth 
curves. In order to apply |[T0] Thm. 5.2.7], we have to study more precisely the behaviour of the operator 
s^d^ + &l at these four vertices. Actually, since the operator e^d^ + &l coincides with its principal part and 
has constant coefficients, we only have to compute the measure wy of the angle at each vertex V of 57 (w). 
Obviously, ijO[±i_-i) = 7r/2 while 

"^"^-°^=""^"n i+o..(±i))0 n"'2 
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Since uoy e (0,7r/2] for V G {(-1, -1), (1, -1), (-1, 0), (1, 0)}, it follows from ^ Thms. 5.2.2 
and 5.2.7] that no singularity occurs at the vertices and that ip e Wp{Q,{u)) for all p G (2,oo). The 
classical Sobolev embedding then implies that t/j G C^~^°' (n(u)) for all a £ (0, 1). Combining this reg- 
ularity with that of u gives that a; i-> \dxtp{x,u{x))\^ + \dztp{x,u{x))\^ belongs to C"([— 1, 1]) and 
Schauder estimates applied to ( 12.5b guarantee that u E C^+"([— 1, 1]) for a G (0, 1). 

Furthermore, since u E C^+"([— 1, 1]), it is easy to checkthat il{u) satisfies an exterior sphere condition 
at each boundary point x e dfl{u). So |]9] Thm. 6.13] appHed to (ll.lOl i- dl.lll i yields ip £ C{i}{u)) n 
C^"*"" (il(ii)). Finally, as T{u) surely is a C^"*"" boundary portion of dn{u), we may invoke [9j Lem. 6.18] 
to deduce that i' e (7^+" (n{u) U r(M)) . □ 



The proof of Theorem ll.ll is thus complete. 



3. The vanishing aspect ratio limit: Proof of Theorem |1.2| 

We now prove Theorem 11.21 and thus consider a family of solutions (we, V'£)ee(o,i) to ( 11.1 Oi l -( fTTST l 
satisfying the bounds (I1.14l i and (11.15b for some fixed A > 0. 
For e e (0, 1), we set 

with from (O and 

<^e{x,v) ■= (t)e{x,v) - V , {x,T])£n. 

We first derive estimates on $e which are uniform with respect to e G (0,1). In the following, K denotes 
an arbitrary positive constant depending only on A and kq. 

Lemma 3.1. There exists a positive constant Ki depending only on A and kq such that, for e G (0, 1), 

imL^in) < 1, (3.1) 

\m\L,in) < K,V^, (3.2) 

Wdri'^eh^m < /Ue, (3.3) 

\\d^,<fe\\L,in) < K,e\ (3.4) 

Proof. Since < c/ij < 1 by (12. 10b . we readily obtain ( 13. lb . We next introduce 
and observe that (11.14b and (11.15b ensure that 

ll/e||L^,(0)< -r + 72 ■ (3.5) 
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Now, it follows from dTHj i- dnjI l that 



/e $e d{x, T]) = 



dxu, 



1 + VI 



d{x,T]) 



1 + eV (9^^,) 



2 / ^a; ''^e 



d(a;,?7) 



2 



1 + Ue y 1 + Me 
2 



$e d(a;,7y) 



(1 + 



+ ^ 2 



^ a?ue 



1 + Me / 1 + Me 



$e 9,,$e d(.T,77) 



We deduce from ( |1.14t . ( 11.151 ), ( 13. It , and the above identity that 



/ /e$ed(x,77) > ||a,,<l>e|| 

Jn 



2 

L2(0) 



Kg Kq / 



|5,,$e||L2(0) 



while ( 13.1b and (13.5b ensure that 



/ /e$,d(.,^)<(^ + £^) 

Jn \ ^0 '^0/ 



Combining the above two inequalities gives ( 13. 3b . Next, thanks to ( 12. 13b and ( 13. lb , we have 

/ |$e(a;,^)|'d(a;,r;) = f f a„$e(a;,y)dy d(x,ry)< / / 
Jn Jn Jri Jn Jo 



d(x,ry)< / / \d,j^,ix,y)\'dydix,7^) 



< |^^|||9„*e"' 



and ( 13.21 ) readily follows from the previous inequality and ( 13.3b . 

Finally, setting Ce and We :— dxdr,^e, we infer from ( I2.12I )-( I2.13I ) that 



1 + £2,^2 (g^y^)2 



in 
Since 



2 ^x'^E 
1+U, 



Ce We 



d(x, 77) = / /e (1 - 5^$e) Ce d(a;, 77) . 

JO 



/ 52$,Ced(x,77)= / d(x,,7) 
Jo Jo 
by GUI Lem. 4.3.1.2 & 4.3.1.3], the above identity also reads 



Ce^ 



(1 + u,y 



+ We - ?7 



1 + Ue 



d(x,77)= / /e (l-9,,$e) Ced(x,77), 



Consequently, owing to ( 11.14b . ( 13.31 ) and ( 13.51) , we deduce from the above identity that 



lie 



< 



< 



Ce^ 



(l + ^^e)^ 

Ke' ||CellL2(n), 



2 d(x,7;)<|i/e|U^(n) (\n\'/' + \\d,'fE\\ 



whence ( 13. 41 ). 



□ 



In order to pass to the limit as e — > in (12.5b , we need to control the trace of 5,,$e on (—1, 1) x {1}. 
For that purpose, the following lemma will be adequate. 
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Lemma 3.2. Given ■& e 14^2^(17), we have 

Proof. We first assume that G C°°(51). For x G (—1, 1) and 7] G (0, 1), we have 

Jri 

Integrating this identity with respect to -q over (0, 1) gives 



d^d{x, 1) 



dri'&{x, rf) dry 



Jii 



d'^d{x,y) dydi]. 



and thus 

\d,A^, 1)1 < \\d,A'X, .)IIl.(o,i) + Wd^^i^, OIU.Co.i) • 
Raising both sides of the above identity to the square and integrating with respect to x over (—1, 1) gives 
the claimed estimate for smooth functions. The general case follows by a density argument, see, e.g., fl] 
Thm. 3.181 or ifTS] Chapt. 2, Thm. 3.11. □ 



A control on the trace of on (—1, 1) x {1} follows at once from Lemma [3T| and Lemma [T2l 
Corollary 3.3. There exists a positive constant K2 depending only on A and kq such that, for e e (0, 1), 

lia^$e(.,i)llL.(-i4) <^2£. (3.6) 

Proof of Theorem \1.2\ By dl.lSI l and the Arzela-Ascoli theorem, there are a sequence (£fc)fc>i and a func- 
tion uq e W'^{-1, 1) such that ^ in (0, 1) and 

— >uo in Wl^{-l,l), 

u,,^uo in Wli^l,!). 

Owing to (11.14b . (13.6b . ( 13.7b . and the definition of $e, we have 

> uo(a;) > Ko- 1, .Te[-l,l], 

and 



1+4 {d^u,,Y 



1 



in Li(-l,l). 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



Combining ( 13. 8b and ( 13. 10b . we may pass to the limit as ej. in ( 12. 5b and conclude that wq is a solution 
to the small aspect ratio equation (11. 9b . where wo(±l) = is guaranteed by ( 12.6b and (13.7b . 
Also, by (ITa . 



and, since 



lim 



> 



1 i-u^^ix) 



1 + z 



1 

1 + z 



dz dx 



1 + 



1 



dz dec , 



we readily obtain (11.161) . where ^jq is given by ( 11.61 ) with u = u^. 
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